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ABSTRACT 

We present new semi-analytical models for the formation of disk galaxies with the purpose of inves- 
Q\ ' tigating the origin of the near-infrared Tully-Fisher (TF) relation. The models assume that disks are 

0^ ! formed by cooling of the baryons inside dark halos with realistic density profiles, and that the baryons 

conserve their specific angular momentum. Adiabatic contraction of the dark halo is taken into account, 
as well as a recipe for bulge formation based on a self-regulating mechanism that ensures disks to be 
O !■ stable. Only gas with densities above the critical density given by Toomre's stability criterion is con- 

sidered eligible for star formation. A Schmidt law is assumed to prescribe the rate at which this gas 
is transformed into stars. The introduction of the star formation threshold density proves an essential 
. ingredient of our models, and yields gas mass fractions that are in excellent agreement with observations. 

Finally, a simple recipe for supernovae feedback is included. We emphasize the importance of extracting 
the proper luminosity and velocity measures from the models, something that has often been ignored in 
the past. We use the zero-point of the if -band TF relation to place stringent constraints on cosmological 
parameters. In particular, we rule out a standard cold dark matter universe, in which disk galaxies are 
too faint to be consistent with observations. The TF zero-point, in combination with nucleosynthesis 
constraints on the baryon density, and with constraints on the normalization of the power spectrum, 
requires a matter density f2o <J 0.3. The observed if -band TF relation has a slope that is steeper than 
simple predictions based on dynamical arguments suggest. Taking the stability related star formation 
threshold densities into account steepens the TF relation and decreases its scatter. However, in order 
[ for the slope to be as steep as observed, further physics are required. We argue that the characteristics 

of the observed near-infrared TF relation do not reflect systematic variations in stellar populations, or 
cosmological initial conditions. In fact, feedback seems an essential ingredient in order to explain the ob- 
served slope of the if-band TF relation. Finally we show that our models provide a natural explanation 
for the small amount of scatter that makes the TF relation useful as a cosmological distance indicator. 

+3 . Subject headings: galaxies: formation — galaxies: fundamental parameters — galaxies: spiral — 

galaxies: kinematics and dynamics — galaxies: structure — dark matter. 

1. INTRODUCTION tions of galaxies. In this paper we use new semi-analytical 

k> ' ,„,,, , , j j c i c , • models for the formation of disk galaxies to investigate the 

Within the standard picture for galaxy formation m Qr Tully-Fisher (hereafter TF) relation (Tully & 

H , which galaxies are thought to form out of the gas which ^. <? 1977) hi h ha th fo 
dissipatively collapses within the potential wells provided 
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by virialized dark halos (White & Rees 1978), disk galaxies 
take a special place. Their flatness and rotational support 

suggest a relatively smooth formation history without vi- Q ^ in tems of a Unear law in a log _ log plane . 
olent non-linear processes. The structure and dynamics 

of disk galaxies are thus expected to be strongly related M = a + &logU ro t, (2) 
to the properties of the dark halos in which they are em- 
bedded. This principle is in fact the key idea behind the where M is the absolute magnitude of the disk, and 
standard model for disk formation which was set out by b = -2.57 is called the slope of the relation. In the light of 
Fall & Efstathiou (1980), and which has since been ex- \,h& strong coupling between disk and dark halo mentioned 
tended upon by numerous investigators (Faber 1982; Fall above, a scaling relation like this is central to theories of 
1983; van der Kruit 1987; Dalcanton, Spergel & Summers galaxy formation. Understanding its origin can put con- 
1997; Mo, Mao & White 1998; van den Bosch 1998). In straints on theories of galaxy formation, and in particular 
this picture, the gas radiates its binding energy but retains n the detailed physics that cause the baryonic component 
its angular momentum, thus settling into a rotationally f a dark halo to be transformed into a luminous galaxy, 
supported disk, the scale length of which is proportional As such, any successful theory of (disk) galaxy formation 
to both the size and the angular momentum of the dark should be able to explain the slope, zero-point, and small 
halo. amount of scatter of this fundamental scaling relation. 

Understanding galaxy formation is intimately linked to Numerous studies in the past have addressed the ori- 

understanding the origin of the fundamental scaling rela- g i n of the TF relation. However, no consensus has been 
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reached, and it is currently still under debate whether the 
origin of the TF relation is mainly governed by initial cos- 
mological conditions (e.g., Eisenstein & Loeb 1996; Avila- 
Reese, Firmani & Hernandez 1998; Firmani & Avila-Reese 
1998a, b), or by the detailed processes governing star for- 
mation (Silk 1997; Heavens & Jiminez 1999) and/or feed- 
back (e.g., Kauffmann, White & Guiderdoni 1993; Cole et 
al. 1994; Elizondo et al. 1999; Natarajan 1999). 

One of the main reasons for this lack of consistency is 
the discordant use of luminosity and rotation measures in 
TF relations. The slope and scatter of the TF relation de- 
pend strongly on the photometric band used to measure 
the luminosities: values of 7 increase from ~ 2.5 in the 
B-band to ~ 4 in the infrared (e.g., Aaronson, Huchra & 
Mould 1979; Visvanathan 1981; Tully, Mould & Aaronson 
1982; Wyse 1982; Pierce & Tully 1988; Gavazzi 1993; Ver- 
heijen 1997). For the rotation velocity a large number of 
different measures have been used: HI line widths from 
single-dish observations, velocities measured from Ha ro- 
tation curves, rotation velocities at a fixed number of disk 
scale lengths, maximum observed rotation velocities, rota- 
tion velocities at the last measured point, and the velocity 
of the flat part of the rotation curve. All these different 
measures yield wildly different TF relations (see Courteau 
1997 and Verheijen 1997 for detailed discussions). There- 
fore it is crucial that one extracts the same quantities from 
the models as the ones used in the empirical relation whose 
origin one seeks to explain. Most previous studies have 
ignored the subtleties associated with these different mea- 
sures, and allowed themselves the "freedom" to compare 
their model predictions with the best fitting TF relation 
available in the literature. 

The other reason why we still lack consensus on the ori- 
gin of the TF relation regards the modeling techniques that 
have been used. Several investigators have used numerical 
simulations (e.g., Evrard, Summers & Davis 1994; Navarro 
& White 1994; Steinmetz & Miiller 1994; Tissera, Lam- 
bas & Abadi 1997; Dommguez-Tenreiro, Tissera & Saiz 
1998; Steinmetz & Navarro 1999; Elizondo et al. 1999). A 
general problem, however, is that the disks that form in 
these simulations have specific angular momenta that are 
more than an order of magnitude lower than those of ob- 
served disks. In addition, most of these simulations ignore 
star formation and feedback, making a direct comparison 
between simulated and observed disks a treacherous enter- 
prise (see Evrard 1997 for a review). Semi-analytical mod- 
eling (SAM) of galaxy formation has been more success- 
ful (White & Frenk 1991; Cole 1991; Kauffmann, White 
& Guiderdoni 1993; Lacey et al. 1993; Cole et al. 1994; 
Heyl et al. 1995; Somerville & Primack 1998). Despite 
reasonable success in reproducing certain TF relations, 
these models are hampered by several shortcomings (see 
discussion in § ||). More recently, Avila-Reese, Firmani & 
Hernandez (1998; hereafter AFH98) and Firmani & Avila- 
Reese (1998a, b; hereafter FA98) presented new SAMs that 
focus only on the formation of disk galaxies. These models 
improve upon previous SAMs in several important ways. 
From a detailed comparison with observed TF relations, 
these authors conclude that the TF relation represents a 
fossil of the primordial density fluctuation field. 

In this paper we present new semi-analytical models for 
the formation of disk galaxies. We improve upon previous 
studies by (i) making a well-motivated choice for the lumi- 



nosity and velocity measures that define the TF relation 
whose origin we seek to understand, (ii) by extracting the 
same measures from the models, and (iii) by using new 
semi-analytical models that focus on the formation of disk 
galaxies and that alleviate several shortcomings of previ- 
ous models. The main aim of this paper is to investigate 
the prevailing mechanism that is responsible for the slope, 
zero-point and small amount of scatter of the TF relation. 
Although our models are, in many ways, complementary 
to those of FA98, we reach a different conclusion. 

2. THE EMPIRICAL TULLY-FISHER RELATION 

One of the most detailed investigations of the empiri- 
cal TF relation is that of Verheijen (1997, hereafter V97), 
who obtained B, R, I, and K band photometry, as well 
as detailed HI synthesis observations of a complete sam- 
ple of (bright) spiral galaxies in the Ursa Major cluster 
(see also Tully et al. 1996). V97 made a detailed com- 
parison of different TF relations using different photomet- 
ric bands and different velocity measures. This study has 
yielded a number of important results: (i) the scatter in 
the TF relation is smallest when using the flat part of the 
HI rotation curve, Vfl a t, as velocity measure, (ii) when us- 
ing the HI line widths, the TF relations become shallower 
and harbor a larger scatter, (iii) the results of Sprayberry 
et al. (1995) and Zwaan et al. (1995), that high surface 
brightness (HSB) and low surface brightness (LSB) spi- 
rals follow the same TF relation, are confirmed, (iv) both 
the scatter and the slope of the TF relation are extremely 
sensitive to selection criteria, and (v) the slope of the TF 
relation changes from b = —6.7 in the £>-band to —10.5 in 
the iC-band. 

Luminosities depend more strongly on extinction, stel- 
lar populations, star formation histories, and metallicities 
going towards bluer passbands. Given the large uncertain- 
ties associated with stellar population models and with a 
proper treatment of extinction by dust, and since the TF 
relation is fundamentally a dynamical law, the most natu- 
ral passbands to use are those in the near-infrared in which 
luminosity is most directly related to the luminous mass 
(see e.g., Gavazzi 1993; Gavazzi, Pierini & Boselli 1996). 
If we further limit ourselves to TF relations that are not 
based on HI line widths, since these are particularly hard 
to interpret in terms of fundamental parameters of the 
galaxy (e.g., Haynes et al. 1999), the available data be- 
comes actually very limited. The only TF relation that is 
based on full near-infrared imaging (rather than aperture 
photometry) and has velocity measures that are derived 
from full HI rotation curves, is the if-band TF relation of 
V97: 

M K - 5 \ogh = 0.683 - 10.5 logU flat , (3) 

Here Vflat is the measured rotation velocity at the flat part 
of the HI rotation curve and h = £/o/(100 kms _1 Mpc _1 ). 
Magnitudes have been converted to Mk — 5 log/i adopting 
a distance to the Ursa Major cluster of 15.5 Mpc, which 
corresponds to Hq = 85kms _1 Mpc _1 (Pierce & Tully 
1988, 1992). In what follows we consider equation (||) the 
'fundamental' TF relation whose origin we seek to under- 
stand. 
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3. MODELING THE FORMATION OF DISK GALAXIES 

In order to investigate whether the origin of the TF 
relation is mainly related to cosmological initial condi- 
tions, or to details regarding the processes associated with 
star formation and/or supernovae feedback, we construct 
new semi-analytical models that focus on the near-infrared 
properties of disk galaxies only. We incorporate a number 
of improvements over previous SAMs, all of which (except 
for FA98) assumed all dark halos to have the same specific 
angular momentum. Since the angular momentum sets, to 
a large degree, the surface density of the disk, and since 
s tar f ormation is a strong function of surface density (see 
§ 3.3), a realistic spread in halo specific angular momenta 
yields a spread in luminosities for galaxies forming in ha- 
los of the same mass. This is especially important when 
investigating for example the scatter in the TF relation. 
Also, we take adiabatic contraction into account (Blumen- 
thal et al. 1986; Flores et al. 1993) and consider the more 
realistic universal halo density profiles, rather than singu- 
lar isothermal spheres, which have been found to poorly 
fit the actual density profiles of dark halos in A-body sim- 
ulations (e.g., Frenk et al. 1988; Efstathiou et al. 1988; 
Dubinski & Carlberg 1991; Navarro, Frenk & White 1995, 
1996, 1997; Tormen, Bouchet & White 1997). We use a 
star formation recipe that is different from previous SAMs, 
and mimic observations of the actual rotation curve of the 
galaxy, rather than to simply assume that the observed 
rotation velocity is equal to the circular velocity at the 
halo's virial radius. As we show below, these latter two 
modifications prove to be of crucial importance. 

The dynamical fragility of disks implies that their mass 
aggregation histories can not have been too violent (e.g., 
Toth & Ostriker 1992). We therefore refrain from calculat- 
ing complete merger histories of the dark halos. The study 
by Navarro et al. (1997) has shown that virialized dark ha- 
los are well described by a universal density profile, whose 
properties can be calculated from the Press-Schechter for- 
malism. In a recent study, Jing (1999) has shown that 
halos with significant amounts of substructure can deviate 
significantly from this universal density profile. However, 
these halos are unlikely to yield disk galaxies, and the uni- 
versal density profiles can thus be used to describe the 
end-state of the mass aggregation histories of halos sur- 
rounding spiral galaxies (see also AFH98). 

As shown by Baugh, Cole & Frenk (1996), the merg- 
ing histories and star formation histories are well coupled. 
Ignoring (small) mergers therefore mainly reflects on the 
details regarding the star formation histories. Since we 
focus on the near-infrared properties, we are insensitive 
to such details, given the fact that mass-to-light ratios in 
the near-infrared depend only very weakly on age and/or 
metallicity (see § [3.5] below) . Furthermore, the quiet mass 
aggregation histories of disk galaxies implies smooth star 
formation histories with little scatter (Baugh et al. 1996). 
Therefore, in the SAMs presented here we do not follow 
the detailed time-resolved formation of disk galaxies, but 
limit ourselves to modeling the time-integrated properties. 

3.1. Properties of the dark matter halos 

Virialized dark halos can be quantified by a mass, Af2oo> 
a radius, r2oo> and an angular momentum, J, originating 
from cosmological torques (e.g., Hoyle 1953). Here r 2 oo is 



defined as the radius inside of which the average density of 
the system is 200 times the critical density of the Universe, 
and M200 is simply the total mass inside r2oo ■ In a series of 
papers Navarro, Frenk & White (1995, 1996, 1997) showed 
that virialized dark halos have universal equilibrium den- 
sity profiles, independent of mass or cosmology, which can 
be well fit by 



p(r) = Pent 



So 



(r/r s )(l + r/r s ) 2 ' 



where 



So = 



200 



3 ln(l + c)-c/(l + c) : 



(4) 



(5) 



with c — r2oo/r s the concentration parameter, and r s a 
scale radius. Given the mass A/200 and the specific cos- 
mology, the concentration parameter c can be calculated 
using the procedure outlined in the Appendix of Navarro 
et al. (1997). In what follows we only consider spherical 
halos, and we refer to halos with the density profile of 
equation (||) as NFW halos. 

The angular momentum of the halo is quantified by the 
dimensionless spin parameter, A, defined by 



A = 



J\E\ 



1/2 



GM. 



5/2 ' 



(6) 



Here E is the halo's total energy. Several studies, both an- 
alytical and numerical, have shown that the distribution 
of spin angular momenta of collapsed dark matter halos is 
well approximated by a log-normal distribution: 



p(A)dA 



1 



<7AV27T 



ln 2 (A/A)\ dA 



(7) 



(e.g., Barnes & Efstathiou 1987; Ryden 1988; Cole & 
Lacey 1996; Warren et al._1992). We consider a distri- 
bution that peaks around A = 0.05 with a\ = 0.6. These 
values are in good agreement with the A-body results of 
Warren et al. (1992), and virtually independent of cosmol- 
ogy (e.g., Lemson & Kauffmann 1999). 

3.2. The formation of disk and bulge 

Once the halo properties are known, we can determine 
the global properties of the disks that form inside them. 
The virialization of the halo heats the baryonic matter to 
the virial temperature. As the baryons start to cool they 
loose their pressure support and settle in a disk. We fol- 
low Fall & Efstathiou (1980) and assume that the specific 
angular momentum of the baryonic material is conserved. 
We further assume that the disk that forms has an expo- 
nential surface density profile. The scale-length of the disk 
is proportional to A and r2oo , and can be calculated using 
equation (28) in Mo et al. (1998). 

Self-gravitating disks tend to be unstable against global 
instabilities such as bar formation. We use the stability 
criterion of Christodoulou, Shlosman & Tohline (1995), 
according to which disks are stable if they obey 



2V C 



(8) 
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with Vd and V c the characteristic circular velocities of the 
disk and the composite disk-bulge-halo system, respec- 
tively 2 . As characteristic velocities we consider the cir- 
cular velocities at R = 3Rd, since Mo et al. (1998) have 
shown that this is the radius where typically the rotation 
curve of the composite system reaches a maximum. We 
set a cr ; t = 0.35 throughout, as suggested for gaseous disks 
(Christodoulou et al. 1995). We have tested that none of 
our results depend significantly on the value of this pa- 
rameter. 

We follow the approach of van den Bosch (1998, 1999) 
and assume that a gaseous disk which is unstable accord- 
ing to criterion (||) transforms part of its disk material into 
a bulge component in a self-regulating fashion, such that 
the final disk is marginally stable (i.e., a = a C rit)- We 
describe the bulge by a sphere with a Hernquist density 
profile: 

Pb{r) = t, 7 — ; 77 > (9) 

Z7r r (r + rb) s 

where rj, is a scale length (Hernquist 1990). Andredakis, 
Pclctier & Balcells (1995) have shown that the effective 
radius r e (defined as the radius encircling half of the pro- 
jected light) is directly related to the total luminosity of 
the bulge by the empirical relation 



M B = -19.75- 2.8 log(r e ) 



(10) 



If we now use Mb — Mr = 4.0, in good agreement with 
the average value for bulges (Peletier & Balcells 1997), and 
r e = 1.8153 rfc (Hernquist 1990), we obtain 



n = 0.71 



10 n L 







0.89 



kpc, 



(11) 



which we use to set the scale length of the bulges in our 
models. Our results do not depend on this simple scaling 
assumption, in practice. The main parameter for the bulge 
is its total mass; changes in its actual density distribution 
are only a second-order effect, and do not influence our 
results. 

3.3. Star formation 

Since the A'-band mass-to-light ratio depends only 
weakl y on the metallicity and age of a stellar population 
(see § 3J3), and since integrated star formation histories 
of disk galaxies are not expected to exhibit strong variety, 
the main parameter of interest is the fraction of the gas 
that is converted into stars over the lifetime of the galaxy. 

In a seminal paper Kennicutt (1989) showed that star 
formation is abruptly suppressed below a critical surface 
density, which is associated with the onset of large-scale 
gravitational instabilities, and which can be modeled by a 
simple Toomre disc stability criterion: 



Scrit(-R) 



k{R) 



3.36 GQ 



(12) 



(Toomre 1964). Here Q is a dimensionless constant near 
unity, (Jgas is the velocity dispersion of the gas, and k is 
the epicycle frequency given by 



k(R) = V2 



V C (R) 



1 



R dVr 



R \~ ' V C (R) dR 

2 Throughout this paper we use the subscripts l d 



This critical density sets the fraction of the gas that is 
eligible for star formation (see also Quirk 1972). Given 
the exponential surface density of the gas disk, with scale 
length Rd and central surface density So, the radius Rsf, 
at which the density of the disk equals S cr j t can be calcu- 
lated. The disk mass inside this radius, and with surface 
density above E C rit is considered eligible for star formation 
(cf. Quirk 1972; Kauffmann 1996). Unless stated other- 
wise we adopt CT gas = 6kms~ 1 and Q = 1.5, as these values 
have been shown by Kennicutt (1989) to yield values of 
i?SF that correspond to the radii at which star formation 
is truncated. 

What remains is to quantify the fraction, /sf, of this el- 
igible mass that is actually turned into stars over the age 
of the galaxy. The star formation rate of disk galaxies is 
well fitted by a simple Schmidt (1959) law: 



V>sfr = a E™ 



(14) 



For a large sample of disk galaxies, Kennicutt (1998) found 
a = 0.25 ± 0.07 Mopc^Gyr -1 and n = 1.4 ± 0.15. If we 
assume that a fraction esf of the mass formed in new stars 
is (instantaneously) returned to the gas phase by means 
of stellar winds and supernovae, one obtains for the gas 
density at the present 



S gas (i?) 



A 



SF 



igal + 



— R/Rd 



1/m 



(15) 



(cf. Heavens & Jiminez 1999), with m = 1/(1 — n), 
Asf = (1 — csf) a, and i ga i the age of the galaxy, which 
we define as the time since the halo's collapse. If we now 
couple this to the threshold density criterion, we can cal- 
culate the total mass in stars formed in the disk. If we 
further assume that all the gas that is transformed into a 
bulge component forms stars at 100 percent efficiency, we 
obtain 



-Rsf 



M* = M b + 2tt 



E c 



-R/Rd 



g{R) 



RdR 



with 



g(R) 



MR) 
rttCR) 



(E g 
(E g 



rit) 
rit) 



(16) 



(17) 



1/2 

■ (13) 

'?>', and 'h' to refer 



Thus, the surface density of the gas is not allowed to drop 
below the critical density given by equation (p2|). 

The total if-band luminosity of the disk-bulge system 
then follows from Lk = M^/T^, where is the mass- 
to-light ratio of the stellar mass. 

3.4. Feedback from supernovae 

The dominant source of feedback that can keep gas out 
of the disk is provided by supernovae (SN). A mass M* of 
stars produces a total amount of energy by SN equal to 

Btotal = Vsn M* £ S n (18) 
to the disk, bulge, and dark matter halo, respectively. 
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Here E'sn is the energy produced by one SN, and 7?sn is 
the number of SN expected per solar mass of stars formed. 
Throughout we adopt t]sn — 4 x 1CP 3 M© -1 , consistent 
with a Scalo (1986) initial mass function (IMF). 

In what follows, we assume that this amount of energy 
can prevent a mass Mhot from ending up in the final disk- 
bulge system, either by means of completely expelling the 
gas from the halo, or by keeping gas in the halo at the virial 
temperature, thus preventing its collapse to the disk. Tak- 
ing account of the escape velocity of the halo, and requiring 
energy balance, one obtains 



Mi 



hot 



M* £sn ??SN Esn 



V 2 

v 200 



(19) 



(cf. Kauffmann et al. 1993; Natarajan 1999), where we 
introduce a new parameter, £sn, which is a measure of the 
efficiency with which the SN energy can prevent the gas 
from ending up in the disk or bulge. One might imagine 
that the free parameter egN is not constant for all galax- 
ies, and we therefore introduce the following simple scaling 
relation 

V200 



250 kms" 1 



(20) 



We can now define the fraction 



M\ 



hot 



= 3.22 eg 



SN' 



( VSTX \ ( Esn_\ ( V200 \ 1 
\ 0.004 Mg- 1 J \^10 51 crg J V 250 kms- 1 J 



(21) 



Thus, in a halo with V200 = 250km s _1 , a mass M* in stars 
can prevent a mass 3.22 £g N M* from settling in the disk 
or bulge. Stellar winds from evolving young stars add to 
the energy released by SN, and we therefore allow values 
of £sn in excess of unity. 

We can no w combine the star formation recipe out- 
lined in § 3.3 with this feedback model to determine the 
masses in the different phases: M*, A/hot, and M co id, 
where the latter is the amount of (cold) gas in the disk 
that is not turned into stars. Mass conservation requires 
M» + Mcoid + M hot = M baryons . Defining e gf as the galaxy 
formation efficiency, which describes what fraction of the 
baryonic mass inside the halo ultimately ends up in the 
disk/bulge, i.e., 



M d + M b = / bar M 2 , 



we can write 



M, 



-1 -1 



hot 



M d + M b 



(22) 



(23) 



This equation can not be solved analytically, as Mi lot , Ma 
and Mb all depend on e g f . However, as we outline in § 3.6 , 
we can solve for e g f using an iterative procedure. 

Since circular velocities depend on whether or not the 
hot gas is expelled from the halo, we introduce the param- 
eter /escj which describes what fraction of A/hot actually 
escapes the halo. The remaining fraction is assumed to 
be at virial temperature and to follow the same density 
profile as the dark matter. It is straightforward to show 
that 

W^V**^ (24) 



This upper limit is reached in the limit where e g f — > 0, 
and thus applies to very low mass systems for which the 
SN feedback is most efficient. For a baryon fraction of 
/bar = 0.05, neglecting / esc thus yields errors on the cir- 
cular velocities that never exceed 3 percent. Lacking a 
physically well-motivated value for / esc and given its in- 
significance for the results presented in this paper, we set 
/esc = 0.5 throughout, i.e., midway its two extremes. 

3.5. Stellar populations 

In order to convert stellar masses to luminosities we need 
to adopt a mass-to-light ratio for the stellar population, 
T* K . In the near- infrared the mass-to- light ratios of stellar 
populations are fairly independent of metallicity and age, 
but these dependencies increase strongly going towards 
bluer passbands (e.g., Maraston 1998a, b). By focusing on 
the near-infrared TF relation we are thus less sensitive to 
uncertainties related to stellar populations (or extinction 
by dust). Using the latest version of the Bruzual & Char- 
lot (1993) population synthesis models, we compute the 
i^-band mass-to-light ratio of a stellar population with a 
Scalo IMF that has been producing stars with a constant 
rate over the past 10 Gyr (a typical value for the age of 
galaxies since collapse). We find a value of T* K — 0.4 and, 
unless stated otherwise, we adopt this value throughout. 

3.6. Constructing a catalogue of model galaxies 

We now combine all the recipes outlined above to con- 
struct a catalogue of model galaxies. After choosing a 
cosmology (see Table [j]) and setting the free parameters in 
our model, we proceed as follows: 

First we randomly assign the halo a mass in the range 
6.3xl0 9 /i _1 M Q < M 200 < 3.6 X10 12 h' 1 M , correspond- 
ing to 30 kms -1 < V200 < 250 kms -1 . We use the pro- 
cedure outlined in the Appendix of Navarro et al. (1997) 
to compute the halo's collapse redshift, z co ii, from which 
we calculate the age t sa \ and the halo concentration pa- 
rameter c. Next we randomly draw a value for the halo's 
spin parameter from the probability distribution of equa- 
tion (@). 

We start by assuming that all the baryons settle in a 
disk (i.e., e g f = 1 and Mb = 0). Taking adiabatic contrac- 
tion into account, we compute the scale-length and central 
surface brightness of the disk (see Mo et al. 1998, and van 
den Bosch 1998 for details). Next we use equation (||) to 
check whether the disk is stable. If a > a C rit, an iter- 
ative procedure is used to compute the bulge mass that 
yields a — a cr it (each iteration step, the procedure for the 
adiabatic contraction is repeated). After convergence, we 
compute M* and Mhot, from which we obtain a new esti- 
mate for the galaxy formation efficiency, e g f. We iterate 
this entire procedure until Mfc+Mcdd+Mhot is equal to the 
total baryonic mass, yielding a completely self-consistent 
disk-bulge-halo model. 

Finally we extract the observational quantities that en- 
ter in the TF relation. The luminosity of the disk in the 
K-b&nd is determined from M* using the mass-to-light 
ratio T* K . As velocity measure, we use the rotation veloc- 
ity at the last measured point of the HI rotation curve of 
the model galaxy. For observed galaxies, this point basi- 
cally coincides with the position at which the HI column 
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density, N[HI], falls below the detectability limit of the ob- 
servations, and is typically of the order of 10 20 cm -2 . Since 
this column density is a projected quantity, the radius of 
the last measured point depends on both the surface den- 
sity of the disk, and its inclination angle with respect to 
the line-of-sight. We therefore draw a random inclination 



angle, 



< i < i r 



and calculate the circular velocity, 



V^bs: at the radius where the projected HI column density 
is equal to 10 20 cm -2 . As mentioned in § ||, the TF relation 
whose origin we seek to understand (equation || ) is based 
on Vflat, i.e., the velocity at the flat part of the rotation 
curve. In Figure [D we show that the velocity at the radius 
where N[HI] = l(r°cm -2 is a good representation of Va a t, 
justifying our choice of velocity measure. 

Throughout we assume that S gas = 1.3Ehi to take ac- 
count of the mass of helium. Unless stated otherwise we 
set ? m in = 45° and i max = 90° as this corresponds to the 
typical values that observers use to select galaxies for their 
TF samples. Since the TF relation of equation (|), which 
we use to constrain our models, is based on spirals of Hub- 
ble type Sb or later, we only allow model galaxies in our 
sample with bulge-to-disk ratios B/D < 0.2. 

4. COSMOLOGICAL CONSTRAINTS FROM THE 
TULLY-FISHER ZERO-POINT 

It is straigthforward to show that simple dynamics pre- 
dict a TF relation with a slope of —7.5 (Dalcanton et 
al. 1997; White 1997; Mo et al. 1998; van den Bosch 1998; 
Syer, Mao & Mo 1999). For the baryon density we adopt 
here (f2bar = 0.0125 h~ 2 ) the predicted TF can be written 
as 



L K = 7.73 x 10 1 



h 3 Q 



0.3 



V* 



200 kms" 



(25) 



(cf. equation [28] in van den Bosch 1998). The observed 
if-band TF relation of equation (|^) in power-law form 
reads 



L K = 5.68 x 10 10 hT 



V fla 



200 kms- 



4.2 



(26) 



Upon equating the zero-points of the observed and pre- 
dicted TF relations at Voon = 200 krns" 1 we obtain 



n h= 1.363 e„ 



0.3 



Vflat 
V20O 



(27) 



This equation immediately suggests that by setting rea- 
sonable limits on e g f, Tk, and Vflat/U2oo, we can place 
strong upper bounds on Qq h. The mass-to-light ratio T k 
in equations ( |25| ) and (^) is defined as the ratio of mass 
over if-band luminosity for the disk. Since not all the 
mass of the disk is transformed into stars, this is not equal 
to the stellar mass-to-light ratio T* K , i.e., 



T K = T 



A' 



Md 



(28) 



and in general Tk > Y* K . Using the Bruzual & Char- 
lot (1993) stellar population models, we find that ;> 

3 Since the predicted and empirical TF relations have different slopes 
depends on the velocity at which we compare the zero-points. For V200 = 
the comparison at V200 = 200 krns -1 we are being conservative 



O.3M /L for a Scalo IMF, at least if the if-band light is 
not dominated by a young (< 1 Gyr) population. The K- 
band mass-to-light ratio depends rather strongly on IMF 
(much more so than on age and/or metallicity) . However, 
most other IMFs that are frequently used, such as the 
Salpeter IMF, have higher mass fractions at the low mass 
end, and therefore yield higher mass-to-light ratios. We 
can thus consider — 0.3 M Q / L Q a conservative lower 
limit. Together with the physical constraint that e g f < 1, 
we obtain 3 

n h < 1.363 f ^ ] ■ (29) 

The ratio Vfl at /V2oo depends on the concentration of the 
dark halo, c, which, for given f2o, fl\ and h, depends on the 
normalization of the power spectrum, erg. Higher values of 
(78 yield more centrally concentrated halos, giving rise to 
higher values of Vfl at /V2oo- Equation ( f29| ) thus yields an 
empirical constraint on a combination of £!o, f^A, h, and 

0-8- 

The constraint of equation (|29|) is over-conservative, in 
that all baryons available in the dark halo are turned into 
stars. However, this leaves no HI gas in the galaxies, con- 
trary to observations. Using the models described in § |[ 
we may derive a more realistic prediction for the TF re- 
lation of equation (|2^). In order to be conservative we 
ignore feedback (i.e., we set e§ N = 0), we assume that all 
the gas above the threshold density for star formation is 
turned into stars with 100 percent efficiency (i.e., Agp is 
such that /sf = 1-0), and we adopt the minimal mass-to- 
light ratio of T^- = 0.3 M / L Q . We thus aim at creating 
maximally bright galaxies. 

In Figure ||| we plot the if-band TF relation for a sam- 
ple of model galaxies in a SCDM cosmology (see Table [[]). 
The parameters of this model, which we refer to as model 
SI, arc listed in Table ||. The empirical TF relation is 
significantly steeper (6 = —10.5 as compared to b — —8.9 
for model SI). At the bright end, the model galaxies are 
too faint by ~ 0.61 mag at V b s = 200 kms -1 . Since 
the parameters of model SI were chosen to ensure max- 
imally bright galaxies, this rules against the SCDM cos- 
mology. For the galaxies in model SI we find an average 
of (Vobs/V2oo) = 1-49 ± 0.11, and from equation ( ^ ) it is 
immediately apparent that the SCDM model is therefore 
ruled out. 

The average concentration of dark halos (and therewith 
the average value of Vobs/U2oo) decreases with decreasing 
o"8- The TF zero-point thus yields an upper limit on ag. 
Results for a variety of (ilo, ft.)-models are presented 
in Figure ||. In a flat universe with J7a = 0, the zero-point 
of the TF relation puts severe limits on the value of erg : for 
h = 0.5 we obtain eg <J 0.08. Such small values are clearly 
inconsistent with the other constraints on erg , and we thus 
conclude that the TF zero-point can be added to a long 
list of observations that rule against a SCDM universe. 
Flat universes with £Iq < 0.3 are far more successful. This 
owes in large to the fact that the universal baryon fraction 
increases with decreasing f^o- 

The values of V o b s /V2Q0 found in our models depend on 
the density profiles of the dark halos, and are thus model 

of —7.5 and —10.5, respectively, the numerical value of the constraint 
: 250 krns" 1 the value of 1.363 reduces to 1.044. Therefore, by making 
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dependent. Although the exact density profiles in the cen- 
tral regions of dark matter halos are still under debate 
(e.g., Moore et al. 1998; Kravtsov et al. 1998), most stud- 
ies now seem to confirm that the outer density profiles fall 
off as r~ 3 . A generic property of such a density profile is 
that Vobs/^200 will be larger than unity, and we therefore 
conclude that the model dependency of our results is only 
weak. Furthermore, if halo densities have central cusps 
that are steeper than for the NFW profiles (as suggested 
by the high resolution simulations of Moore et al. 1998), 
our estimates of Vobs/V2oo are too low, and our limits on 
us conservative. 

4.1. Comparison with previous work 

The result that the SCDM cosmology yields galaxies 
that are too faint for their rotation velocities has already 
been pointed out by numerous previous studies. However, 
the analysis presented here is very much different, and 
yields important new results. 

Already the very first applications of semi-analytical 
modeling revealed problems with the normalization of the 
TF relation (Cole 1991; White & Frenk 1991; Lacey & 
Silk 1991; Lacey et al. 1993; Kauffmann et al. 1993; Cole 
et al. 1994; Heyl et al. 1995). If the star formation and 
feedback parameters of these models were tuned to fit the 
luminosity function of galaxies, the zero-point of the TF 
relation turned out too faint by ~ 2 orders of magnitude. 
This was interpreted as indicating that there are too many 
dark halos of a given velocity. 

Here we have shown that if one adheres to the nucleosyn- 
thesis constraints on the baryon density, it is in principle 
impossible to fit the TF zero-point for a SCDM cosmology 
with a reasonable value for trg: even after setting the model 
parameters to yield maximally bright galaxies, they turn 
out too faint, independent of whether or not these models 
fit the galaxy luminosity function. The constraints pre- 
sented here are thus more strict than those from previous 
SAM studies. 

Contrary to all previous SAMs, Somerville & Primack 
(1998) were able to simultaneously reproduce the normal- 
ization of the TF relation and the galaxy luminosity func- 
tion in a SCDM cosmology (with a% = 0.67). The reason 
for this apparent contradiction with our results is two- 
fold. First of all, Somerville & Primack used a baryon 
density that is 1.5 times the value adopted here. Sec- 
ondly, they considered isothermal halos, ignored adiabatic 
contraction, and made no distinction between V200 and 
the actual velocity measure that enters the empirical TF 
relation (the same is true for all other previous SAMs). 
We have shown that a proper treatment of the different 
velocities is a crucial aspect of any attempt to try and 
comprehend (the origin of) the TF relation. In particular, 
from equation (|2^) it is immediately apparent that had we 
directly associated V200 with the rotation measure V h s we 
would not have been able to rule against SCDM. 

In addition to the SAMs, several studies based on nu- 
merical simulations have encountered similar problems 
with the normalization of the TF relation, e.g., Steinmetz 
& Navarro (1999) and Elizondo et al. (1999). These au- 
thors suggest that the angular momentum catastrophe and 
declining star formation rates, respectively, are likely to 
be the cause for this discrepancy. The analysis presented 
here, however, shows that the problem is more fundamen- 



tal than that. 

5. CONSTRAINTS FROM THE SLOPE OF THE 
TULLY-FISHER RELATION 

Simple dynamical arguments predict a TF relation with 
a slope of —7.5, very different from the observed value of 
— 10.5. From equation (^5|) it is immediately clear that 
one can "tilt" the TF relation to the observed slope by 
satisfying the condition 

ft (^)V»r (30, 

Each of the three parameters on the left side of this equa- 
tion represents a different class of physics: the galaxy for- 
mation efficiency, e g f, is determined mainly by processes 
related to feedback, the mass-to-light ratio T k is related 
to details regarding the star-formation and stellar pop- 
ulations, and V^oo/Kibs is related to the density of the 
dark halo, which in turn reflects cosmological initial con- 
ditions. Here we investigate which of these processes is 
most important in tilting the TF relation towards its ob- 
served slope. Based on the discussion in the previous §, 
we limit ourselves here to our fiducial ACDM3 cosmology 
(see Table |). 

5.1. Star formation 

The mass-to-light ratio T k in equation (|30| ) is given by 
equation (p8|). The ratio M^/M* is set by the parameters 
Q, which sets the threshold density for star formation, and 
by n and Agp , which control the star formation efficiency. 
In Figure [| we plot the TF relations for three models that 
only differ in the value of Q (see Table ||) . All three models 
have fsF = 1.0 and ignore feedback (i.e., e g f = 1). If we set 
Q = 15.0 (i.e., ten times its fiducial value), S C rit becomes 
so low that virtually all the disk material is converted into 
stars, independent of V^oo- Since V2oo/K>bs depends only 
very weakly on V200 and e g f = 1 we thus obtain a TF re- 
lation with a slope of —7.5. We do not consider a value 
of Q = 15.0 physical, but merely show the results of this 
model for comparison. Setting Q to its fiducial value of 1.5, 
T k is found to decrease with increasing mass. This results 
in a significantly steeper TF relation with b = —9.1, and 
suggests that a further decrease of Q might actually yield 
a TF relation with b ~ —10.5. Indeed, for Q = —0.8 we 
obtain b = —10.0, consistent with the observations within 
the error-bars. Besides influencing the slope of the TF re- 
lation, Q also affects the amount of scatter. For values of 
Q <; 0.8 we can obtain even steeper TF relations, but at 
the cost of introducing unacceptable levels of scatter. The 
origin of this scatter is discussed in detail in § ^. 

We also examined the influence of changing the star for- 
mation parameters Agp and n (such as to obtain fgp < 1) 
while setting Q to its fiducial value of 1.5. In each case, 
we find that the Schmidt law causes /sf to decrease with 
increasing mass, resulting in a TF relation that is shal- 
lower than for model LI. In addition, reducing /sf tends 
to increase the scatter. 

5.2. Stellar populations 

Next we experiment with a mass dependence of the 
stellar mass-to-light ratio T^. Rather than adopting 
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T* K = 0.4 M / L , we set oc (V2oo) s an d determine 
the value of s that yields b = —10.5, while setting all other 
parameters to that of model LI. The best fitting model 
(L3) is plotted in the left panel of Figure [| and has 



r* K 



0.35 



250 kms- 



-0.7 



(31) 



This implies if-band mass-to-light ratios for stellar pop- 
ulations in systems with V200 = 50 kms -1 to be a fac- 
tor three larger than for systems with V200 = 250 kms -1 . 
Stellar population models suggest that such a large spread 
in T* K cannot be due to variations in ages and/or metal- 
licities, but it might reflect systematic variations of the 
IMF with galaxy mass (see e.g., Maraston 1998a). How- 
ever, there is no observational support for such a varia- 
tion. Furthermore, if low mass systems have considerably 
higher if-band mass-to-light ratios than the more mas- 
sive systems, this effect would be much stronger in bluer 
passbands. One would thus predict a TF relation in the 
£>-band with b <C —10.5, in clear contradiction with ob- 
servations. We therefore consider a systematic variation 
of T* K with mass an unlikely explanation for the steepness 
of the observed if-band TF relation. 

5.3. Cosmological initial conditions 

The ratio V2oo/K>bs m equation ( J30| ) depends on the 
concentration c of the dark halo, which is set by the halo's 
collapse redshift z co i\. For a CDM power spectrum, the 
mass dependencies of z C oii, and henceforth c, are found to 
be fairly small (e.g., Navarro et al. 1997). Here we decou- 
ple the calculation of z co \\ from the CDM power spectrum 
and adopt the simple parameterization 



1 + Zcoii = 3.5 



Mo, 



10 12 Mr. 



-c 



(32) 



with £ a free parameter. The constant 3.5 is chosen to fit 
the TF zero-point for halos with M 200 = 10 12 M Q . 

We find that for ( ~ 0.26 the slope of the model TF 
relation becomes as steep as observed. This model (L4) is 
plotted in the middle panel of Figure |[ From a comparison 
of the mass dependence of c for this model with that of a 
number of scale free power spectra in Navarro et al. (1997), 
we infer that £ = 0.26 corresponds roughly to a power 
spectrum of density fluctuations with an effective slope at 
the scale of galaxies of n e s — —1.0. This is significantly 
shallower than for CDM (for which n e g lies between —2.0 
and -2.4 for systems with 10 9 M < Ai 20 o < 10 12 M ), 
and therefore implies a cosmology with more power on 
small scales as compared to CDM. This, however, seems 
inconsistent with observations which imply that, if any- 
thing, the power on small scales has to be less than for 
CDM (see e.g., Kauffmann et al. 1993; Cole et al. 1994; 
Klypin et al. 1999; Moore et al. 1999). 

We therefore conclude that the slope of the observed TF 
relation does not reflect cosmological initial conditions, in 
clear contrast to the results of AFH98. This inconsistency 
is a reflection of the different choices for the empirical TF 
relation used to constrain the models: if AFH98 had cho- 
sen to use the TF relation of V97 as empirical constraint, 
they also would have concluded that the slope of the TF re- 
lation does not reflect cosmological initial conditions. The 



empirical relations used by AFH98 have slopes of b ~ —8 
and are based on HI line widths and if-band aperture pho- 
tometry. Subtle effects can cause non-linearities between 
aperture luminosities and total luminosities, which is not 
taken into account in the modeling of AFH98. In addi- 
tion, HI lincwidths are difficult to interpret in terms of 
fundamental parameters of the galaxy (see § ||). The if- 
band TF relation used to constrain our models is based 
on total luminosities derived from CCD imaging, and on 
rotation measures derived from the complete HI rotation 
curves. We extract the same measures from our model 
galaxies ensuring a fair comparison between models and 
observations. We therefore feel confident that we can rule 
out initial cosmological conditions as the primary cause for 
the observed slope of the TF relation (at least for a CDM 
power spectrum). 

5.4. Supernova feedback 

Having ruled out stellar populations and cosmological 
initial conditions, we finally resort to SN feedback to try 
and obtain a TF relation as steep as observed. 

The right panel of Figure || plots the TF relation of 
model L5, for which we have tuned the feedback parame- 
ters by fitting the empirical if-band TF relation (see Ta- 
ble H for the parameters). The model fits the empirical 
TF relation of V97 extremely well. We experimented with 
a variety of models with different values of £g N and v and 
found that other combinations yield virtually equally good 
fits, whereby lower values of £g N require a steeper mass de- 
pendence (i.e., smaller v). In general, we find that for the 
ACDM3 cosmology and = 0.4M Q /L Q the percentage 
of baryons prevented from ending up in either the disk or 
bulge has to increase from ~ 10 percent for systems like 
the Milky Way (i.e., V ohs = 220 kms" 1 ) to ~ 70 percent 
for low mass systems with V b s = 50 kms -1 . The par- 
ticular feedback model used here induces a curvature in 
the TF relation (cf. Natarajan 1999), but for v ^ -2 
this has a negligible effect over the magnitude range of 
-19 > M K - 51ogft > -25 (see § |). Cole et al. (1994), 
on the other hand, used v = —3.5, resulting in a strongly 
non-linear TF relation. 

5.5. Different cosmological models 

The discussion above is based on the ACDM3 cosmol- 
ogy. Although this model is consistent with the normal- 
izations from both COBE and the cluster abundances, it 
predicts halo concentrations with c ~ 9 on the scale of 
galaxies. Such high values, however, yield rotation curves 
that are too steep to be consistent with several observed 
rotation curves (Navarro 1998). We therefore now consider 
the ACDM2 cosmology (see Table |l), which was shown by 
Navarro (1998) to be consistent with all rotation curves in 
his sample. Halos in this cosmology typically have c ~ 3. 

The panels on the left in Figure || show the results for a 
ACDM2 cosmology with the same parameters as for model 
LI (see Table ||). Galaxies are typically ~ 2 mag too bright 
for their rotation velocities. This owes to both the high 
baryon fraction in this cosmology (/bar = 0.25) and to the 
fact that halos are less concentrated such that Vobs/^oo is 
smaller (i.e., we find (Vobs/V^oo) = 1.26±0.09 as compared 
to (V ohs /V 20Q ) = 1.47 ± 0.11 for model LI). In the mid- 
dle panel, feedback is added with its parameters tuned to 
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fit the TF relation. Within our simplistic feedback model 
we require ~ 200 percent of all the SN energy to prevent 
on the order of 80 to 95 percent of the baryons from set- 
tling in the disk/bulge system. It seems unlikely that the 
combined effect of SNe and stellar winds can be so effi- 
cient (i.e., see Mac Low & Ferrara 1999). We can lower 
the feedback efficiencies while still fitting the TF relation 
by lowering the star formation efficiency Agp (such that 
fsF < !)■ In fact, we find a considerable amount of free- 
dom in the parameters describing SN feedback and star 
formation to obtain good fits to the empirical TF relation. 
In general, however, one requires either a higher feedback 
efficiency, or a lower star formation efficiency than for the 
ACDM3 cosmology. 

6. CONSTRAINTS FROM THE SCATTER OF THE 
TULLY-FISHER RELATION 

The models in the previous sections contain two sources 
of scatter in the TF relations: the spread in halo spin pa- 
rameters, ax, and the spread in inclination angles, which 
we quantify by Ai = i max - i min . 

For a halo of a given mass, a non-zero o\ yields a spread 
in both the magnitudes and rotation velocities of disks. 
Lower values of A yield more compact disks with a larger 
fraction of their mass eligible for star formation (i.e., galax- 
ies are relatively bright). In addition, the compactness of 
the disk induces, because of the adiabatic contraction, a 
more strongly concentrated halo, and therefore a relatively 
high value of V ohs /V 20 o- Model L0 (for which Q = 15.0) 
in Figure ^ depicts a large amount of scatter. Its origin 
is illustrated in the small inset in Figure ^, which shows 
the "TF relation" for model galaxies with the same mass, 
but different spin parameters. Because of the high value of 
Q, the critical density for star formation is relatively low. 
This implies that the radius i?sF where star formation is 
truncated occurs at relatively large radii of ~ 8i?d, where 
the growing curve of the cumulative disk mass is very flat. 
Consequently, scatter in A translates mainly into scatter 
in Vobs and not in Mk ■ The opposite applies to model L2, 
for which Q = 0.8. Here, Rsf ~ 3Rd, and a\ translates 
mainly into scatter in disk magnitudes. Most remarkably, 
for our fiducial value of Q (i.e., model LI), o\ translates 
into scatter in both Mk and Vobs but such that galaxies 
with different values of A are aligned with the actual TF 
relation. Consequently, o\ contributes only weakly to the 
scatter in the TF relation. For Q = 1.5 we thus obtain 
a natural explanation for the small amount of scatter ob- 
served. 

Different values of i for otherwise identical galaxies re- 
sults in different radii where N[HI] = 10 20 cm~ 2 , and thus 
in different values of Vobs- As shown in Figure |, the ro- 
tation measures used in the model TF relations are fairly 
insensitive to i, indicating that the rotation curves are 
close to flat. Consequently, Ai only adds a small fraction 
of the total scatter. 

An additional source of scatter, which has hitherto not 
been taken into account, stems from the variation in mass 
aggregation histories (MAHs), which induces scatter in the 
density profiles of halos of the same mass. For halos that 
are well represented by a NFW density profile, this im- 
plies a scatter in the concentration parameter c. Thus if 
we know how scatter in MAHs translates into scatter in c 
we can investigate how this cosmological scatter adds to 



that of the TF relation. This information is provided by 
the high resolution numerical simulations of Jing (1999). 
In our models c is computed from the collapse redshift of 
the halo, and we can thus introduce scatter in c by insti- 
tuting scatter in z co ii- To that extent we introduce a new 
parameter, a z , which sets the variance in collapse redshifts 
of a given mass, and tune it to reproduce the scatter in c 
in the simulations of Jing. Although the amount of scatter 
increases strongly with the amount of halo substructure, 
the fragility of disks suggests that they are embedded in 
relatively smooth halos. We therefore normalize a z by fit- 
ting the distribution of concentration parameters in halos 
without significant amounts of substructure. In Figure |^ 
we plot the log-normal distribution that best fits the nu- 
merical simulations of Jing, together with a histogram for 
halos in model L5 and with a z — 0.4. Except for a small 
offset in the zero-point, which is not important for ob- 
taining an appropriate measure for the amount of scatter, 
the agreement is very good. We can thus use a z — 0.4 
as a representative measure for the standard deviation in 
MAHs for halos surrounding disk galaxies. This translates 
into scatter in the concentration of halos, and thus in the 
ratios Vobs/^200 and the surface densities of the disks. 

In order to investigate the main source of scatter in the 
TF relation, we compare four TF relations that only differ 
in the values of Ai, o\ and a z (see Figure ||). The mod- 
els have the same parameters as model L5. The upper 
right panel, for which Ai = a\ = a z = clearly reveals 
the small amount of curvature eluded to above, and which 
causes a non zero scatter of <jm — 0.06 mag around the 
best fitting linear relation. The main source of this curva- 
ture is the SN feedback. Setting a\ to its fiducial value of 
0.6 increases the scatter to <7m = 0.13 mag. If we further 
add a variation in inclination angles, with Ai = 45°, the 
scatter becomes that of model L5: cfm — 0.16 mag. If 
finally we set a z to its fiducial value of 0.4, this further 
increases gm to 0.20 mag. We find that for o z ;> 0.4, 
the dispersion in halo collapse redshifts becomes the main 
source of scatter in the TF relation. 

V97 derived confidence levels on the amount of intrin- 
sic scatter in the TF relation of equation (^), and found 
cm < 0.30 mag at 95 percent confidence level, with a most 
likely value of 0.11 mag. This is in excellent agreement 
with our value ctm = 0.20 mag for a z = 0.4. 

7. SCATTER FROM COSMOLOGICAL INITIAL CONDITIONS 

Eisenstein & Loeb (1996; hereafter EL96) used the 
excursion-set formalism of Lacey & Cole (1993) to com- 
pute MAHs for halos of a given mass, and calculated pre- 
dictions for the amount of scatter in the TF relation result- 
ing from this variation in MAHs. They concluded that the 
amount of scatter in the observed TF relations is smaller 
than what one expects from the scatter in MAHs. How- 
ever, we have shown in § [| that our models predict an 
amount of scatter that is in good agreement with observa- 
tions, even after we take the variance in MAHs, modeled 
quantitatively by cr z , into account. Here we investigate the 
reason for this disagreement. 

In their analysis, EL96 assume that the mass gained in 
an infinitesimal redshift interval is accreted in a smooth 
and spherical fashion, and calculate the induced change in 
binding energy, E. Using the assumption that halos are 
isothermal, E cx V200, and any variation in MAHs thus 
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produces scatter in the TF relation. EL96 assume that 
the velocity measure used in the TF relation is strictly de- 
termined by V200 and that the overall luminosity of the 
galaxy is directly related to the mass of its halo. For a 
SCDM model with ct§ = 0.65 they find om — 0.36 mag, 
and conclude that this is too large to be consistent with 
the observed scatter. 

We can make a rough comparison between our method 
and that of EL96, using the fact that <r z can be used to pa- 
rameterize scatter in MAHs. The analysis of EL96 is based 
on the spherical collapse model which, for halos with the 
same mass and in an Einstein-de Sitter universe, yields 
E oc (1 + z co \\). Combining this with the assumption of 
isothermal halos we thus obtain V b s oc (1 + ^coii) 1 ^ 2 - In 
our analysis, a z induces scatter in halo concentration pa- 
rameters, which in turn yields variance in V^bs/V^oo- The 
calculation of c from z co \\ is based on the Press-Schechter 
formalism, and is well calibrated against high resolution 
A-body simulations. In Figure [| we compare the dis- 
tributions of V hs that results from scatter in MAHs de- 
rived using both methods. The method of EL96 predicts 
a standard deviation in rotation velocities that is approx- 
imately a factor two larger than for our method. If we 
halve the amount of scatter expected by EL96, we obtain 
om = 0.18, bringing their results in good agreement with 
ours. Since our method does not make any assumptions 
regarding spherical collapse, uses more realistic halo den- 
sity profiles, and has been well calibrated, we feel confident 
that our estimates of the expected scatter in the TF re- 
lation are the more reliable. We therefore disagree with 
EL96, and conclude that the observed amount of scatter 
in TF relations is consistent with expectations from cos- 
mological initial conditions. 

8. CONSTRAINTS FROM GAS MASS FRACTIONS 

Our results regarding the TF relation would be of little 
use if the model galaxies were not representative of real 
galaxies. Here we briefly examine the gas mass fractions 
of our model galaxies, as they provide useful, additional 
constraints. 

McGaugh & de Blok (1997; hereafter MB97) have ob- 
tained gas mass fractions, defined by 



M, 



cold 



M cold + M* 



(33) 



for a sample of 108 spiral galaxies. In Figure [h] we plot 
/gas as function of absolute magnitude and central sur- 
face brightness for the data of MB97 and for three of 
our models that all fit the empirical TF relation of V97. 
Model L5 agrees extremely well with the observations, in- 
dicating that the observed gas mass fractions in spirals 
are consistent with being governed by a star formation 
threshold density that originates from Toomre's stability 
criterion with Q = 1.5. Model L2, for which Q = 0.8, 
yields gas mass fractions that are slightly too large. To- 
gether with the fact that this model predicts a TF zero- 
point that depends on surface brightness (see § [}|) and 
that Q — 0.8 is inconsistent with the empirically deter- 
mined value prompts us to rule against this model. Model 
L8, for which /sf < 1, predicts gas mass fractions that 
are clearly too large to be consistent with observations. 
We find this to be a general problem for models in which 



only a fraction of the mass eligible for star formation is 
actually turned into stars over the lifetime of the galaxy. 
Virtually all models with Q and A$f set to their fiducial 
values (that have been determined empirically) yield gas 
mass fractions in good agreement with the data of MB97, 
suggesting that indeed fsF — 1-0. 

9. SUMMARY AND CONCLUSIONS 

We have used new semi-analytical models to investigate 
the origin of the near-infrared TF relation. By focusing 
on the near-infrared we are less susceptible to uncertain- 
ties related to extinction by dust and stellar populations. 
As empirical constraint we have used the A-band TF re- 
lation of V97, which is the only near-infrared TF relation 
based on CCD imaging (rather than aperture photome- 
try) combined with full HI rotation curves (rather than 
HI line widths). An important improvement over previous 
studies aimed at understanding the TF relation concerns 
our special care in extracting the same luminosity and ro- 
tation measures from our model, as the ones used in the 
empirical TF relation. 

Since dark halos have density profiles that at large radii 
fall off more rapidly than an isothermal (i.e., the NFW 
profiles used here fall off as r~ 3 ), the circular velocities of 
the disk at the flat part of the rotation curve are gener- 
ally larger than the circular velocities at the virial radius. 
This has important consequences for the TF zero-point, as 
it allows us to put constraints on cosmological parameters. 
In a SCDM universe, dark halos are too centrally concen- 
trated, and the baryon fraction is too low, to be able to fit 
the A-band TF relation. The observed zero-point puts an 
upper limit on the normalization of the power spectrum 
of as <, 0.08, in clear contradiction with constraints from 
COBE and the abundances of rich clusters of galaxies. 
We have shown that the observed TF zero-point (com- 
bined with the nucleosynthesis constraints on the baryon 
density) favor a Universe with f2o ^ 0-3- Previous studies 
of the TF relation based on SAMs were only able to rule 
against SCDM, if the models were normalized to fit the lu- 
minosity function of galaxies. The reason why we can rule 
out SCDM without this additional constraint is our more 
sophisticated treatment of the relationship between V200 
and V ohs . 

Whereas simple dynamics predict a TF relation with a 
slope of —7.5, the A'-band TF relation reveals a slope of 
— 10.5. This implies that the physics regulating star for- 
mation and feedback, coupled with the mass dependence 
of halo densities and stellar populations, has to tilt the 
TF relation to its observed slope. The introduction of a 
stability-related star formation threshold density increases 
the slope of the TF relation, reduces its scatter, and yields 
gas mass fractions that are in excellent agreement with 
observations. Setting Toomre's Q parameter to its empir- 
ically determined value tilts the TF relation to a slope of 
~ — 9, and additional physics are thus required to obtain 
a TF relation as steep as observed. We have presented 
four different physical mechanisms that all yield good fits 
to the observed A-band TF relation: lowering Q to ~ 0.8, 
systematic variations in the A-band mass-to-light ratios, 
a power spectrum of initial density fluctuations that dif- 
fers from that for CDM, and feedback. Except for the 
latter, each of these results in inconsistencies with other 
observations, and can be ruled out as the prevailing mech- 
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anism for the observed characteristics of the near-infrared 
TF relation. 

The feedback efficiency required to tilt the TF relation 
to its observed slope has to be higher in lower mass sys- 
tems, and depends strongly on Qq and h (if the baryon 
fraction is taken to be constrained by the nucleosynthesis 
results). In particular, in a flat universe with f2o = 0.2 
and h — 0.5 (a cosmological model that yields low density 
halos in agreement with observed rotation curves), ~ 200 
percent of the available SN energy is required to fit the 
observed TF zero-point. Although in principle unphysi- 
cal, we do not wish to draw too strong conclusions from 
our oversimplified model for feedback. We do emphasize, 
however, that some amount of feedback is required to yield 
TF relations with a slope as steep as observed. More so- 
phisticated models of SN feedback, such as in the work 
of Mac Low & Ferrara (1999), are required to investigate 
whether our inferred feedback efficiencies are realistic. Our 
best fitting models with feedback reveal a small amount of 
curvature in the TF relation, in qualitative agreement with 
recent observations of dwarf galaxies (Mathews, van Dricl 
& Gallagher 1998; Stil & Israel 1998). Further constraints 
on the models come from the gas mass fractions, which 
require that virtually all the gas mass eligible for star for- 
mation (i.e., with densities above the critical density) is 
transformed into stars over the lifetime of the galaxy. 

For the fiducial value of Q the scatter in halo spin pa- 
rameters affects both the luminosities and the observed 
rotation velocities, but to such a degree that galaxies are 
scattered along the TF relation, rather than perpendicular 
to it. The introduction of star formation threshold densi- 
ties thus yield a natural explanation for the small amount 
of scatter observed. Taking account of a realistic spread 
in halo concentrations, consistent with what one expects 
from the spread in mass aggregation histories, our model 
that best fits the zero-point and slope of the if -band TF 
relation predicts a scatter of <~ 0.2 mag only, in excellent 
agreement with observations. In a follow-up paper (van 
den Bosch & Dalcanton 1999) we compare the models pre- 
sented here to numerous other observational constraints. 

This work has benefited greatly from discussions with 
Marcella Longhetti, Claudia Maraston, Jeroen Stil, and 
Marc Verheijen. I am indebted to Julianne Dalcanton and 
George Lake for advice and critical assessments of an ear- 
lier draft of the paper, to Stephane Chariot for provid- 
ing results from his stellar population models, and to the 
anonymous referee for his comments that helped improve 
the presentation of the paper. Support for this work was 
provided by NASA through Hubble Fellowship grant # 
HF-01102.11-97.A awarded by the Space Telescope Science 
Institute, which is operated by AURA for NASA under 
contract NAS 5-26555. 
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Fig. 1. — The ratios V o ^ s /V200 (left panel) 
with the parameters of model L5 (see Table 



ad R bs/Rd (right panel) as functions of the cosine of the inclination angle for model galaxies 
I. Here i? Q b s is defined as the radius where the projected surface density of the HI reaches a 
column density of 10' !u cm — 2 , and V b s = V c (R bs)- For disks seen close to edge-on [i ~ 90°), i? bs becomes as large as ~ 10Rd- For close 
to face-on systems, R h a /Rd — 4.5 with a considerable scatter. Despite this relatively large change in -R bsi the ratio V o b s /V200 is almost 
independent of i. This indicates that Vob s is an accurate measure of the flat part of the rotation curve, and we therefore use V b s and Vfl at 
without distinction throughout the paper. 




log[V J (km/s) 

Fig. 2. — Model TF relation for model SI (SCDM with parameters set to obtain maximally bright galaxies). Open circles correspond to 
the model galaxies, and the solid line corresponds to the observed K-hand TF relation of V97. Only galaxies with —19 > Mk — 51og?i > —25 
are shown, which corresponds to the magnitude interval on which the empirical TF relation is based. The slope of the model TF relation is 
too shallow (b = —8.9 as compared to 6 oba = —10.5), and model galaxies are too faint at the bright end. Since galaxies are maximally bright, 
this rules against the SCDM model (see text). 
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Fig. 3. — Constraints on erg and h for four flat cosmologies with different values of Qo- The light gray areas correspond to the parts of 
parameter space allowed by the TF zero-point. Cosmological models with (crg,h) in the white areas yield galaxies that are too faint. In 
addition, constraints are plotted from the 4 year data from the COBE DMR experiment (Bunn & White 1997; solid lines) and from the 
cluster abundances at z = 0: ct 8 = (0.52 ± 0.04) n-°-52+0.l3« ( Ek6j Cole & Frenk 1996; dark gray bars ) and z ~ 7: CTg = ( g ± .1) Q~ - 29 
(Bahcall & Fan 1998; median gray bars). In a flat universe without cosmological constant (upper right panel) the TF zero-point imposes 
stringent constraints on erg and h. In order for galaxies in this cosmology to be as bright as required by the TF relation, <rg < 0.08 when 
h = 0.5. This is in clear contradiction with the other constraints. Since a decrease of Qo induces an increase of the baryon fraction (i.e., we 
adopt / bar = 0.0125 Q^h^ 2 throughout) , the TF constraints become less restrictive for lower values of the matter density. Combining all 
constraints suggests Qq < 0.3. 
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i.B 3 2.2 2.4 1.8 2 3,2 3,4 1,8 3 3.2 2.4 

l°e[V g J (km/s) 

Fig. — The upper panels show TF relations for three models in a flat ACDM universe with Qq = 0.3 (symbols are the same as in 
Figure El). All these models ignore feedback (i.e., e g f = 1) and have all the mass that is eligible for star formation turned into stars over the 
lifetime of the galaxy (i.e., /sf = 1). The models only differ in the value of Q as indicated in the panels. In addition, the panels list the model 
ID and the slope of the best fitting TF relation. The small insets plot TF relations for model galaxies that only differ in A (all galaxies have 
M200 = 3 X lO 11 /! -1 Mq and have a value of A drawn from the log-normal distribution of equa tion |7| with o\ = 0.6). These insets show how 
a\ contributes to the scatter (see discussion in § p|). The lower panels plot log(T/f ) (equation |28|) as function of log(V t, s ). The dotted lines 
correspond to Tjj- = = 0.4 Mq/ Lq. Lower mass systems require a significantly higher value of Tjf than more massive galaxies in order 
to fit the observed slope. This can be accomplished by lowering Q to ~ 0.8. 




Fig. 5. — Model TF relations for three ACDM3 models that all fit the empirical A-band TF relation of equation (||). Each of these models, 
however, represent a different class of physics. In model L3, the A-band stellar mass-to-light ratios are varied with mass to obtain b = —10.5. 
In model L4, the same is accomplished, but by adopting a power spectrum that is shallower than for CDM on the scale of galaxies. Finally, 
in model L5 feedback is included to tilt the TF relation to its observed slope. 
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Fig. 6. — Model TF relations for three ACDM2 models (upper panels). The models differ only in the parameters describing star formation 
and feedback (see Table G). The effects of these two physical processes are indicated in the panels in the two lower rows, where we plot e g f 
(a measure of the feedback efficiency) and /gp (a measure of the star formation efficiency) as functions of log(V t, s ). ff no feedback is taken 
into account, and all the mass eligible for star formation is turned into stars over the lifetime of the galaxy (panels on the left), galaxies are 
almost two magnitudes too bright. Models L7 and L8, in which feedback is taken into account, both provide good fits to the observed TF 
relation. 
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0.5 1 1.5 

c/c 

Fig. 7. — Histogram of the normalized distribution of c/cq for model L5 with cr z = 0.4. Here cq is the halo concentration parameter 
predicted by the Prcss-Schechter formalism (sec Navarro ct al. 1997 for details). The solid curve represents the log- normal distribution of 
c/co for halos without significant amounts of substructure as determined by Jing (1999) from high resolution iV-body simulations (in Jing's 
notation, we only included halos with dvi max < 0.15). Despite a small offset in the mean (Jing finds (c) ~ 0.95co), the two distributions 
agree well, and we conclude therefore that a standard deviation of 0.4 in halo collapse redshifts is a good representation of the scatter in halo 
concentrations that owes to scatter in mass aggregation histories. 




L5 S 2.5 3 3,5 

log[V e J (km/a) log[V g J (km/s) 

Fig. 8. — Illustration of the build-up of scatter in the TF relation of model L5. Contrary to previous plots (i.e., Figures 0, Q El and pj) 
we plot the TF relation over the entire range of Mjf of our model galaxies. This emphasizes the curvature present in the TF relation, and 
which owes mainly to the feedback. In the upper left panel, no intrinsic scatter is present. The small amount of scatter in the TF relation 
of om = 0.06 mag owes to the curvature. The upper right panel shows the same TF relation but now with a\ = 0.6. The variance in halo 
spin parameters increases ctm to 0.13 mag. Setting Ai = 45° further increases the scatter to om = 0.16 mag (lower left panel). Finally, the 
addition of scatter in the halo collapse redshifts (<r z = 0.4) yields a total scatter of ctm = 0.20 mag. 
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0.6 0.8 1 1.2 1.4 



Fig. 9. — The shaded histogram plots the distribution of "observed" velocities (normalized to the mean) for model SI, but with Ai = ax = 
and a z = 0.5. The thick solid line corresponds to the same distribution, but calculated from the procedure used by Eisenstein &: Loeb (1996). 
This latter method is based on the simplified assumptions of spherical accretion and isothermal halos, and predicts roughly twice the amount 
of scatter in the TF relation. 



Frank C. van den Bosch 



19 




-14 -16 -18 -3D -22 26 S4 S3 20 
M B -51ogh TO ^ (mag/arc sec E ) 



Fig. 10. — The gas mass fraction / ga a (see equation |33| ) as function of absolute B-band magnitude (left panels) and central disk surface 
brightness (right panels). Upper panels plot the data oi McGaugh & de Blok (1997), while the panels in the lower three rows correspond to 
models L2, L5 and L8, as indicated. The data is based on S-band photometry, and we convert our models (K-ba,nd) to the i?-band using 
B - K = 0.64 - 0.18(Mjf - 51ogh 75 ) with h 75 = H /(75 kms -1 Mpc -1 ) (see Appendix A in van den Bosch & Dalcanton 1999). The thin 
solid lines have no physical meaning, but are plotted to facilitate a comparison between models and data. Model L5 yields gas mass fractions 
that are in excellent agreement with the data. The model reproduces the decrease of / ga s with increasing luminosity and surface brightness, 
and the absolute values of /g aB are in excellent agreement with the data of MB97. Models L2 and L8, on the other hand, predict gas mass 
fractions that are too large to be consistent with observations. 
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Table 1 

Parameters of cosmological models. 



Model 


n 


n A 


h 


0-8 


/bar 


to 


Clusters 


COBE 


RCs 


TF 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


SCDM 


1.0 


0.0 


0.5 


0.6 


0.05 


13.0 


Y 


N 


N 


N 


ACDM3 


0.3 


0.7 


0.7 


1.0 


0.085 


13.5 


Y 


Y 


N 


Y 


ACDM2 


0.2 


0.8 


0.5 


0.38 


0.25 


21.0 


N 


Y 


Y 


Y 



Note. — Column (1) lists the name of the CDM cosmology by which we refer to it in the text. Columns (2), (3), 
(4), and (5) list the matter density, f^o, the density in the form of a cosmological constant, S1a = A-/{3Hq) (with A the 
cosmological constant), the hubble parameter, h = 7?o/(100kms _1 Mpc -1 ), and the rms mass variance on a scale of Shr 1 
Mpc, as- Column (6) lists the baryon fraction /bar = ^bar/^o based on a baryon density of f2bar = 0.0125 hr" 1 as required 
by nucleosynthesis constraints (e.g., Walker et al. 1991). Column (7) lists the age of the universe in Gyr. Columns (8) 
to (11) indicate whether the models are consistent with the observed cluster abundances, the COBE normalization, the 
observed rotation curves of spirals (RCs; see Navarro 1998), and the TF zero-point, respectively (see § for details) 



Table 2 

Parameters of models discussed in the text. 



ID. 


Cosmology 


Q 




n 


A SF 


£ SN 


p 


C 


b 


cm 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0) 


(10) 


(11) 


SI 


SCDM 


1.5 


0.3 


1.4 


0.25 


0.0 


0.0 




-8.9 


0.15 


L0 


ACDM3 


15.0 


0.4 


1.4 


0.25 


0.0 


0.0 




-7.5 


0.28 


LI 


ACDM3 


1.5 


0.4 


1.4 


0.25 


0.0 


0.0 




-9.1 


0.18 


L2 


ACDM3 


0.8 


0.4 


1.4 


0.25 


0.0 


0.0 




-10.0 


0.23 


L3 


ACDM3 


1.5 


(x (V 2 oo)- ' 7 


1.4 


0.25 


0.0 


0.0 




-10.5 


0.20 


L4 


ACDM3 


1.5 


0.4 


1.4 


0.25 


0.0 


0.0 


0.26 


-10.5 


0.22 


L5 


ACDM3 


1.5 


0.4 


1.4 


0.25 


0.05 


-0.3 




-10.5 


0.16 


LG 


ACDM2 


1.5 


0.4 


1.4 


0.25 


0.0 


0.0 




-9.0 


0.19 


L7 


ACDM2 


1.5 


0.4 


1.4 


0.25 


2.0 


0.0 




-10.6 


0.13 


L8 


ACDM2 


1.5 


0.4 


1.4 


0.01 


0.05 


-2.6 




-10.7 


0.13 



Note. — Column (1) lists the model ID, by which we refer to it in the text. Column (2) lists the cosmological model 
used (see Table |TJ). Columns (3) and (4) list the value of the Toomre parameter, Q, and the if-band stellar mass-to- 
light ratio, T* K , respectively. The star formation parameters n and Ag-p are listed in columns (5) and (6), respectively. 
Columns (7) and (8) list the feedback parameters £g N and v, respectively. Column (9) gives the parameter £. This 
parameter is only defined for model L3; for all other models, the mass dependence of c is determined from the CDM 
power spectrum. Finally, columns (10) and (11) list the slope b and scatter ctm (in mag) of the TF relation for model 
galaxies with — 19 > Mk — 5\ogh > —24 and with a z = 



